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Abstract 

The possibility of construction of continuous analogue of Poisson 
distribution with the search of bounds of confidence intervals for parameter 
of Poisson distribution is discussed. Also, in the article is shown that the 
true value of a parameter of Poisson distribution for the observed value x 
has Gamma distribution with the scale parameter, which is equal to one, 
and the shape parameter, which is equal to x. The results of numerical 
construction of confidence intervals are presented. 
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1 Introduction 

In article the unified approach to the construction of confidence intervals 
and confidence limits for a signal with a background presence, in particular 
for Poisson distributions, is proposed. The method is widely used for the 
presentation of physical results |2| though a number of investigators criticize 
this approach (in particular, this approach avoids a violation of the cov- 
erage principle). Series of Workshops on Confidence Limits has been held 
in CERN and Fermilab. At these meetings demands for properties of con- 
structed confidence intervals and confidence limits have been formulated [^. 
On the other hand, the results of experiments often give noninteger values of 
a number of observed events (for example, after background subtraction [^) 
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when a Poisson distribution take place. That is why there is a necessity 
to search a continuous analogue of Poisson distribution. The present work 
offers some generalization of Poisson distribution for continuous case. The 
generalization given here allows to construct confidence intervals and con- 
fidence limits for Poisson distribution parameter both for integer and real 
values of a number of observed events, using conventional methods. More 
than, the supposition about continuous of some function /(x,A) described 
below allows to use a Gamma distribution for construction of confidence 
intervals and confidence limits of Poisson distribution parameter. In present 
article we consider only the construction of confidence intervals. 

In the Section 2 the generalization of Poisson distribution for the contin- 
uous case is introduced. An example of confidence intervals construction for 
the parameter of analogue of Poisson distribution is given in the Section 3. 
In the Section 4 the results of construction of confidence intervals having the 
minimal length for the parameter of Poisson distribution using a Gamma 
distribution are discussed. The main results of the article are formulated in 
the Conclusion. 



2 The Generalization of Discrete Poisson Distri- 
bution for the Continuous Case 

Let us have a random value ^, taking values from the set of numbers x € X. 
Let us consider two-dimensional function 

f{x,X) = ^e-\ (1) 

where x > A > 0. Assume, that set X includes only integer numbers, then 
discrete function f{x, A) describes distribution of probabilities for Poisson 
distribution with the parameter A and a random variable x 

Let us rewrite the density of Gamma distribution using unconventional 
notation 

where a is a scale parameter, x > — 1 is a shape parameter]^ and A > is a 
random variable. Here the quantities of x and A take values from the set of 

conventional notation for the shape parameter is M r = a; + 1, r > 0. 
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real numbers. Let us set a = 1 and denote xl = r(x + 1), then a continuous 
function 



f{x,X) = —e'\ X>0, x>-l (3) 

is the density of Gamma distribution with the scale parameter a = 1. 

To anticipate a little, it is indicative of the Gamma distribution of pa- 
rameter A for the Poisson distribution in case of observed value x = x. 

Figure 1 shows the surface described by the function /(x,A). Smooth 
behaviour of this function along x and A (see Fig. 2) allows to assume that 
there is such a function Z(A) > — 1, that 

f{x,X)dx = l (4) 

for given value of A. It means that in this way we introduce continued ana- 

A^ 

logue of Poisson distribution with the probability density f(x,X) = — rc"'** 

xl 

over the area of the function definition, i.e. for x > /(A) and A > 0. The 
values of the function f(x, A) for integer x coincide with corresponding mag- 
nitudes in the probabilities distribution of discrete Poisson distribution. De- 
pendences of the values of function /(A), the means and the variances for 
the suggested distribution on A have been calculated by using programme 
DGQUAD from the library CERNLIB 0] and the results are presented in 
Table 1. This Table shows that series of properties of Poisson distribution 
{E^ = A, = A) take place only if the value of the parameter A > 3. 
It is appropriate at this point to say that 

fi^, A)dx = / dx = e-V(A). (5) 

7o r(x + l) 



The function 







is well known and, according to ref. ||8|, 

m= E wf-Yy + 0(lA|-^-''-^) = e^ + 0(|A|-^) (7) 

n=~N ^ ' 
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if A — oo, |ar(7A| < ^ for any integer N. Nevertheless we have to use the 
function l(X) in our calculations in Section 3. We consider it as a mathe- 
matical trick to illustrate a possibility of construction of confidence intervals 
by numerically in the case of real value x. 

In principle, we can numerically transforms the function /(x, A) in the 
interval x S (0, 1) so that 

/■oo roo 

/ /(x, X)dx = 1, E^= xf{x, \)dx = A, 
Jo Jo 

poo 

Di= / {x-Eiff{x,\)dx = \ 
Jo 

for any A. In this case we can construct confidence intervals without intro- 
ducing of /(A). 

In Section 3 assumption about continuous of the function /(x,A) along 
the variable x are used for construction of confidence intervals of a parameter 
A for any observed x. 

Let us construct a central confidence intervals for the continued analogue 
of Poisson distribution using function /(A). 



3 The Construction of the Confidence Intervals 
for Continued Analogue of Poisson Distribution. 

Assume that in the experiment with the fixed integral luminosity (i.e. a 
process under study may be considered as a homogeneous process during 
given time) the x events of some Poisson process were observed. It means 
that we have an experimental estimation A(x) of the parameter A of Pois- 
son distribution. We have to construct a confidence interval (Ai(x), A2(x)), 
covering the true value of the parameter A of the distribution under study 
with confidence level 1 — a, where a is a significance level. It is known from 
the theory of statistics [^], that the mean value of a sample of data is an 
unbiassed estimation of mean of distribution under study. In our case the 
sample consists of one observation x. For the discrete Poisson distribution 
the mean coincides with the estimation of parameter value, i.e. X = x. This 
is not true for a small value of A in the considered case (see Table 1). That 
is why in order to find the estimation of A(x) for small value x there is nec- 
essary to introduce correction in accordance with Table 1. Let us construct 
the central confidence intervals using conventional method assuming that 

|_"/(x,Ai)dx = | (8) 
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for the lower bound Ai and 




for the upper bound A2 of confidence interval. 

Figure 3 shows the introduced in the Section 2 distributions with param- 
eters defined by the bounds of confidence interval (Ai = 1.638, A2 = 8.493) 
for the case x = \ = A and the Gamma distribution with parameters a = 1, 
X = X = A. The association between the confidence interval and the Gamma 
distribution is seen from this Figure. The bounds of confidence interval with 
90% confidence level for the parameter of continued analogue of Poisson dis- 
tribution for different observed values x (first column) were calculated and 
are given in the second column of the Tabic 2. It is necessary to notice that 
the confidence level of the constructed confidence intervals always coincides 
exactly with the required confidence level. As it results from Table 2 that 
the suggested approach allows to construct confidence intervals for any real 
and integer values of the observed number of events in the case of the val- 
ues of parameter A > 3. The Table 2 shows that the left bound of central 
confidence intervals is not equal to zero for small x, what is not suitable. 

Also note that 90% of the area of Gamma distributions with parameter 
X = X are contained inside the constructed 90% confidence intervals for 
observed value x. However, for small values of x we have got values of the 
area close to 88%, i.e. less than 90%. 



4 Shortest Confidence Intervals for Parameter of 
Poisson Distribution. 

As is follow from the formula (3) for f{x,X) (sec Fig. 3) the true value of 
A of Poisson distribution for the observed value x has Gamma distribution 
(at least both for integer x > if we consider pure Poisson distribution 
and for real .x > if we consider analogue of Poisson distribution) with the 
parameters a = 1 and x = x, i.e. 

P{X\x) = ^e-\ (10) 
x] 

This supposition allows to choose confidence interval of a minimum length 
from all possible confidence intervals of given confidence level without vio- 
lation of the coverage principle. The bounds of minimum length area, con- 
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taining 90% of the corresponding Gamma distribution square, were found 
numerically both for integer value of x and for real value of x. We took 
into account that x = x and found the shortest 90% confidence interval for 
the parameter of Poisson distribution. The results are presented in third 
column of the Table 2. For comparison with the results of conventional 
procedure Q of finding confidence intervals, the results of calculations of 
confidence intervals for integer value of x Q are adduced in the Table 2. 
By this means confidence intervals, got using Gamma distribution, may be 
used for real values of x. 

5 Conclusion 

In the article the attempt of introducing of continued analogue of Poisson 
distribution for the construction of classical confidence intervals for the pa- 
rameter A of Poisson distribution is described. Also, in the article is shown 
that the true value of a parameter of Poisson distribution for the observed 
value X has Gamma distribution with the scale parameter a = 1 and the 
shape parameter x = x. Two approaches (with using of function and 
with using of Gamma distribution) are considered. Confidence intervals for 
different integer and real values of number of observed events for Poisson 
process in the experiment with given integral luminosity are constructed. 
The second approach allows to construct any confidence intervals for a pa- 
rameter of Poisson distribution. 
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Table 1: The function 1{X), mean and variance versus A. 



A 


/(A) 


moan (-B^) 


variance (-D^) 


0.001 


-0.297 


-0.138 


0.024 


0.002 


-0.314 


-0.137 


0.029 


0.005 


-0.340 


-0.130 


0.040 


0.010 


-0.363 


-0.120 


0.052 


0.020 


-0.388 


-0.100 


0.071 


0.050 


-0.427 


-0.051 


0.113 


0.100 


-0.461 


0.018 


0.170 


0.200 


-0.498 


0.142 


0.272 


0.300 


-0.522 


0.256 


0.369 


0.400 


-0.539 


0.365 


0.464 


0.500 


-0.553 


0.472 


0.559 


0.600 


-0.564 


0.577 


0.653 


0.700 


-0.574 


0.681 


0.748 


0.800 


-0.582 


0.785 


0.844 


0.900 


-0.590 


0.887 


0.939 


1.00 


-0.597 


0.989 


1.035 


1.50 


-0.622 


1.495 


1.521 


2.00 


-0.639 


1.998 


2.012 


2.50 


-0.650 


2.499 


2.506 


3.00 


-0.656 


3.000 


3.003 


3.50 


-0.656 


3.500 


3.501 


4.00 


-0.647 


4.000 


3.999 


4.50 


-0.628 


4.500 


4.498 


5.00 


-0.593 


5.000 


4.997 


5.50 


-0.539 


5.500 


5.497 


6.00 


-0.466 


6.000 


5.996 


6.50 


-0.373 


6.500 


6.495 


7.00 


-0.262 


7.000 


6.995 


7.50 


-0.135 


7.500 


7.494 


8.00 


0.000 


8.000 


7.993 


8.50 


0.000 


8.500 


8.496 


9.00 


0.000 


9.000 


8.997 


9.50 


0.000 


9.500 


9.498 


10.0 


0.000 


10.00 


9.999 
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Table 2: 90% C.L. intervals for the Poisson signal mean A for total events 
observed x. 





bounds 


(Section 3) 


bounds 


(Section 4) 


bounds 


(ref[l]) 


X 


Ai 


A2 


Ai 


A2 


Ai 


A2 


0.000 


0.121E-08 


2.052 


0.0 


2.303 


0.00 


2.44 


0.001 


0.205E-08 


2.054 


0.0 


2.304 






0.002 


0.292E-08 


2.056 


0.0 


2.306 






0.005 


0.666E-08 


2.061 


0.0 


2.311 






0.01 


0.307E-07 


2.076 


0.0 


2.320 






0.02 


0.218E-06 


2.098 


0.0 


2.337 






0.05 


0.765E-05 


2.166 


1.66E-05 


2.389 






0.10 


0.137E-03 


2.275 


2.23E-05 


2.474 






0.20 


0.186E-02 


2.490 


6.65E-05 


2.642 






0.30 


0.696E-02 


2.692 


1.49E-04 


2.806 






0.40 


0.161E-01 


2.891 


2.60E-03 


2.969 






0.50 


0.295E-01 


3.084 


5.44E-03 


3.129 






0.60 


0.466E-01 


3.269 


1.35E-02 


3.290 






0.70 


0.673E-01 


3.450 


2.63E-02 


3.452 






0.80 


0.911E-01 


3.629 


4.04E-02 


3.611 






0.90 


0.1179 


3.804 


6.12E-02 


3.773 






1.0 


0.1473 


3.977 


8.49E-02 


3.933 


0.11 


4.36 


1.5 


0.3257 


4.800 


0.2391 


4.718 






2.0 


0.5429 


5.582 


0.4410 


5.479 


0.53 


5.91 


2.5 


0.7896 


6.340 


0.6760 


6.220 






3.0 


1.056 


7.076 


0.9284 


6.937 


1.10 


7.42 


3.5 


1.340 


7.792 


1.219 


7.660 






4.0 


1.638 


8.493 


1.511 


8.358 


1.47 


8.60 


4.5 


1.946 


9.188 


1.820 


9.050 






5.0 


2.264 


9.869 


2.120 


9.714 


1.84 


9.99 


5.5 


2.590 


10.55 


2.453 


10.39 






6.0 


2.924 


11.21 


2.775 


11.05 


2.21 


11.47 


6.5 


3.264 


11.87 


3.126 


11.72 






7.0 


3.609 


12.53 


3.473 


12.38 


3.56 


12.53 


7.5 


3.961 


13.18 


3.808 


13.01 






8.0 


4.316 


13.82 


4.160 


13.65 


3.96 


13.99 


8.5 


4.677 


14.46 


4.532 


14.30 






9.0 


5.041 


15.10 


4.905 


14.95 


4.36 


15.30 


9.5 


5.406 


15.73 


5.252 


15.56 






10. 


5.779 


16.36 


5.640 


16.21 


5.50 


16.50 


20. 


13.65 


28.49 


13.50 


28.33 


13.55 


28.52 
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Figure 1: The behaviour of the function f{x,X} versus A and x if f{x,X) < 1. 
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Figure 2: Two-dimensional representation of the function f{x,X) versus A and x 
for values f{x,X) < 1. 
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Figure 3: The probability densities f{x, A) of continuous analogue Poisson distri- 
bution for A's determined by the confidence limits Ai and A2 in case of observed 
number of events x — -i and the probability density of Gamma distribution with 
parameters a = 1 and x ~ x ~ 4. 
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